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Abstract 

We study integrated and differential rates for the production of charged Higgs 
bosons of the Minimal Supersymmetric Standard Model via 6-quark initiated 
subprocesses in pp collisions at the Large Hadron Collider. In detail, we compute 
cross sections and distributions of the reactions: bU — > bDH + — > bDr + v T ® C.C. 
and bU — > bDH + — > bDtb — > bDbb}] © C.C, for a scalar in the intermediate 
(i.e., < mt + m,},) and heavy (i.e., M ± > mt+m^) mass range, respectively (U 
and D represent generic u- and d-type light quarks). In the first case, a detailed 
treatment of various possible backgrounds is also given. 
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1. Introduction 



At the Large Hadron Collider (LHC) |l], 0, the charged Higgs boson H ± of the Minimal 
Supersymmetric Standard Model (M.SSM.) is expected (if it exists) to be copiously 
produced in top quarks decays, via the chain t — > bH ± — > b(ru T ), provided that m t > 
Mh± + rrib and that the value of tan f3 is low or high enough^- 

The main decay modes of the H scalar are (for M H ± < m t + rrib) into cs, tv t 
and W ± h pairs, where in the last case h — ► bb (h being the lightest neutral Higgs 
boson of the M.SSM.) j3]. The first decay is never dominant, whereas the second is 
overcome by the third in a narrow mass window right before the opening of the tb Higgs 
decay threshold, but only at very small values of tan (3. Otherwise, the branching ratio 
BR(H ± — > rv T ) is the largest (around 98%, for tan/5 > 2), and it depends only slightly 
on the f3 angle. When Mjj± > m t + m,b, the H — > tb decay mode is in practise the 
only relevant one (with a branching ratio of practically 100%)Q 

Extensive studies and simulations for H ± production at the LHC have been carried 
out 0, ||. Top quarks are produced in ti pairs, via qq and gg fusion, with a large cross 
section (around 500 pb at \fs = 14 TeV and for m t = 175 GeV). The signal that has 
been considered in the ATLAS and CMS Technical Proposals is the one involving one 
top decaying to a charged Higgs, and the other decaying inclusively (i.e., either via a 
or, mostly, a W^) into electrons and muons (and corresponding neutrinos). The 
charged Higgs boson is searched for by means of the leptonic signature — > tv t . 

Since neutrinos in the final state prevent one from reconstructing the Higgs mass 
from the momenta of its decay products, the existence of H ± signals in the data can 
be inferred only from an excess of r production with respect to what is predicted for 
the Standard Model (SJA) backgrounds (lepton universality breaking signal) . Among 
the letter, one must number the irreducible ones (non-resonant r production and ti 
production followed by t — W^b and W ± — > tv t ) as well as the reducible ones (mainly 
ti where either a jet from a W ± fakes a r or a 6-jet decays leptonically into or fakes a 
r, but also bb production followed by bb — > r + jets and W ± + jets, with one of the jets 
faking a r) [p], 0. By selecting an isolated high p T lepton and by requiring one jet to 
have a high transverse energy together with one 6-tagging should allow one to explore 
a large portion of the (Ma, tan f3) plane, with a significance up to 5a (assuming a 10 
fb -1 integrated luminosity of the collider) |]J. 

It is the purpose of this letter to study other production mechanisms of charged 
Higgs bosons of the M.SSM. at the LHC, via subprocesses with 6-quarks in the initial 
state. In particular, we calculate (for a H ± scalar whose mass is below the tb threshold) 
the signal reaction (including charge conjugation) 

bU -> bDH + -> bDr + v T ® C. C, (1) 

and the background processes 

bU -> bDr + v T © C. C, via intermediate W ± , 7, Z, H, h, A, (2) 

1 The minimum of the t — > bH^ decay rate is at about tan/3 = 6. 

2 In the above discussion and throughout this paper we have assumed that the mass scale of the 
Supersymmetric particles is above the mass, such that only decays into ordinary matter are here 
considered. 



bU -> bDr + u T © C.C., via intermediate W ± ,g, (3) 

[/£) -> 66r + z/ r © C.C., via intermediate W^g, (4) 

where represents a generic u(d)-type light quark (i.e., u, d, s and c) found inside 

the proton and H, h, A are the three neutral Higgs bosons of the AASSM.. For a heavy 
we consider the production and decay chain (again including charge conjugation) 

bU -> 6M+ -> -> 6Z>66jj © C. C, (5) 

assuming that, because of the spectacular signature that is produced in the final state, 
background processes can easily be kept under control. In reactions (0)-(3) and (0) we 
treat the initial 6-quark as a constituent of the proton with the appropriate momen- 
tum fraction distribution fb/ P (x,Q 2 ), as given by our default set of partonic structure 
functions. 

The relevance of these reactions can be understood if one considers that: 

• at the typical (partonic) energies of the LHC the content of 6-quarks inside the 
colliding protons is very much enhanced, compared to lower energy hadronic 
scatterings (such as, e.g., at the Tevatron); 

• the presence of 6-quark in initial, virtual and final states means that many of the 
M.SSAA Yukawa couplings of the Higgs bosons of the theory are increased for 
large values of tan j3; 

• the vertex tagging performances of the LHC detectors are expected to become 
almost 'ideal' by the time the machine starts operations, thus allowing one to 
greatly reduce the QCD background of light quark and gluon jets. 

Finally, we also stress that 6-quark initiated processes at LHC energies have already 
been demonstrated to be important in the case of neutral AiSSM. Higgs production, 
especially at large tan/3's and for intermediate masses of the scalars M, as well as in 
the case of charged Higgs production, via the reaction gb — > H + t © C.C. ||. 

The plan of this paper is as follows. In the next Section we give some details of the 
calculation and list the values adopted for the various parameters. Section 3 is devoted 
to a discussion of the results. Conclusions are in Section 4. In the Appendix we write 
down the amplitudes squared of the signal processes. 



2. Calculation 

To calculate processes (0)-© we have used the spinor techniques described in Refs. || 
0, The FORTRAN codes we have produced have been counter-checked against the 



outputs of MadGraph 0, which incorporates the HELAS subroutines flO |. We have 
always found perfect agreement between the two kind of programs. The codes written 
using the helicity formalism of Refs. || ^, U have also been tested for gauge invariance. 
Furthermore, in order to speed up the numerical evaluations in Monte Carlo simulations, 
the matrix elements for the signal processes (0) and (H) have been computed using the 
textbook method of taking the trace of the gamma matrices, with the help of FORM 



flTTf . The analytical expressions obtained in this way are very simple, so that we do 
reproduce them here (see the Appendix). They have been eventually implemented and 
their numerical results agree with those of the other codes. The integrations over the 



phase space have been performed using VEGAS [[[4 . 

The tree-level Feynman diagrams that one needs for computing processes (|]]) - @ 
are given in Fig. la,b,c and dQ. The labelling of the particles in the figures corresponds 
to their ordering in the left-hand side (for the initial state) and in the right-hand side 
(for the final state) of equations dJ)-©. Note that the virtual particle content of the 
diagrams is explicitly indicated in Fig. 1 for all reactions. 

Concerning the values of the various parameters entering in the computation of 
processes (|l])-(||), we have proceeded as follows. First, we have set up the mass scale 
of the Supersymmetric partners of ordinary matter well above the energy reach of the 
LHC, such that we can neglect their contribution in our calculations. To further simplify 
the discussion, we have assumed a universal soft Supersymmetry-breaking mass [JT^, |TJ]] 

m l = m l = m % ( 6 ) 
and negligible mixing in the stop and sbottom mass matrices, 

A t = A b = /1 = 0. (7) 

One-loop corrections to the masses of the M.SSM neutral CP-even Higgs bosons 
and to the mixing angle a are introduced via the parameter e of Ref. [TE\, given by 
(neglecting the 6-mass) 

rnfl»(l + ^, (8) 



" 8tt 2 M^ ± sin 2 6 W 
where e 2 = 47ra em . One then gets JT3 



± {[(Ml - M 2 Z ) cos 2f3 + e/ sin 2 f3} 2 + [M\ + Mf ) 2 sin 2 2/?} V2 , (9) 

(M\ + M 2 )sin2/3 

For the M.SSM. charged Higgs masses we have maintained the tree-level relations 

M^ ± =M 2 +M^ ± , (11) 

since one-loop corrections are small compared to those for the neutral Higgses jyj . 

In the numerical calculations presented in the next Section we have adopted the 
following values for the electromagnetic coupling constant and the weak mixing angle: 
a em = 1/128 and sin 2 6w = 0.2320. The strong coupling constant a s , which appears at 
next-to-leading order in the computation of the charged Higgs decay width (see Ref. 0) 

3 The C.C. diagrams can be obtained by simple crossing and time inversion of fermion lines. 



and 



and enters in some of the production mechanisms, has been evaluated at two loops, with 
= 230 MeV, and with the number of active flavours Nj (and the corresponding 

A^|^) calculated according to the prescription of Ref. [[16] at the scale Q 2 = s. 

For the gauge boson masses and widths we have taken Mz = 91.1888 GeV, Y z = 2.5 
GeV, My/± = 80.23 GeV and 1TV± = 2.08 GeV, while for the fermion masses we have 
used m e = m Ue = = m VT = 0, = 0.105 GeV, m T = 1.78 GeV, m u = = 
m s = m c = 0, rrif, = 4.25 GeV and m t = 175 GeV, with all widths equal to zero except 
for T t . We have calculated this at tree-level within the AiSSAi, using the expressions 
given in Refs. [[17], The universal Supersymmetry-breaking squark mass is in the 
numerical analysis rriq = 1 TeV and the LHC centre-of-mass (CM) energy is \^s pp = 14 
TeV. Finally, throughout the paper we have always used MRSA pOf as the default set 
of partonic distributions, with the same a s and A^- as above. 

3. Results 

As it is impractical to cover all possible regions of the AiSSAi parameter space 
(Ma, tan (3), we have decided to concentrate here on the two representative (and ex- 
treme) values tan/? = 1.5 and 30, and on masses of the pseudoscalar Higgs boson A in 
the range 60 GeV ~ Ma ~ 500 GeV. The large bibliography existing on the AiSSAA 
Higgs decay phenomenology allows one to easily extrapolate our results to other values 
of tan (3 frf|. 

In Fig. 2 we display the total cross sections at the LHC for processes (l)-(4), for 
values of Ma up to 480 GeV and for the two above-mentioned tan/3's (note that a 
minimum transverse momentum of 10 GeV is required for all detectable particles in 
the final states of all processes). The relevant feature of Fig. 2 is that process (1) has 
rather large rates, of the same order as the backgrounds (2)-(4). This is particularly 
true below Ma ~ 120 GeV (which corresponds to M^± « 145 GeV, small window in 
the upper right corner of Fig. 2), and more for large than for small values of tan (3. The 
latter aspect is a consequence of the fact that the BR of the charged Higgs boson into 
tv t pairs is enhanced at tan j3 = 30, a value for which the cs channel is negligible (see 
Ref. ||). An additional contribution comes from graph 2 in Fig. la, which involves 
Yukawa vertex contributions proportional to tan (3. However, the largest part of the 
signal cross section is due to graph 1 in Fig. la, because of a resonant top decay. The 
steep decrease of the signal rates around Ma ~ 150 GeV is due to the opening of the 
H ± — > tb off-shell decay channel (see Refs. |3], |22||). 

Concerning the background processes (2)-(4), one notices that they are 'roughly' 
independent of Ma and tan/?. This is obvious for processes (3)-(4), as they proceed 
through SAi graphs (Figs, lc and d), whereas for process (2) this indicates that the 

4 Actually, we have done done so only in the production processes (those represented by the graphs 
in Fig. la,b,c and d). In fact, for process (5), to describe the decay chain t — > bW^ — > 6jj, we have used 
a Narrow Width Approximation (NWA) for the top, by implementing the decay formulae as given in 
Ref. . This has been done in order to avoid a large consume of CPU time in computing exactly a 
2 — > 6 partonic process convoluted with initial structure functions and in presence of multiple resonant 
peaks in different regions of the phase space. We are confident that such an approximation does not 
spoil the validity of our conclusions. 



contributions to the total cross section due to interactions involving M.SSM. vertices 
(that is, graphs 8 and 11 in Fig. lb) are irrelevant (even at tan /3 = 30). We also notice 
that, unlike processes (2)-(3), reaction (4) is not an irreducible background, as its final 
state is different from that of the signal. Nevertheless, as it includes two bottom quarks 
among the produced particles, the probability that its final signature be the same as 
the signal is (neglecting the misidentification of light flavour quarks as fe's as well as 
correlations between the two possible tags): P = 1 — B, where B = 1e h — e\, being 
the efficiency of tagging one 6-jet. Hence, better 6-tagging leads to lower detection rates 
of the background process (4). 

Therefore, from the figure it is clear that in principle a large excess of r events could 
be produced in the scattering process, provided that Ma - 120 — 130 GeV (at large 
tan/3's such an interval can be possibly extended up to 130-140 GeV). The possibilities 
of actually disentangling the signal depend strongly on the detector performances of the 
LHC, in particular in recognising the displaced vertex in jets which originate from b- 
quarks^. In fact, the signature that one would look for is 6jrF/ T X, where i^r represents 
the missing (transverse) energy due to the neutrino escaping the detectors and j is 
the jet arising from the light parton scattered in the proton. In order to quantify the 
significance of the signals we list in Tab. I the total cross sections of processes (1)- 
(4) (as read from Fig. 2) for M A = 60(80) [100]{120} GeV (corresponding to M H ± = 
100(113) [128]{144} GeV), for both values of tan/5 = 1.5 and 30, multiplied by the 
6-tagging efficiencies and by the yearly luminosity / Cdt = 10 fb _1 (the minimum 
considered for the final collider design). Furthermore, we also compute the relative 
excess of signal events S respect to the total background B, as A = S/B x 100. For 
the microvertex performances, we have adopted the following three reference values: 
Sb = 1.0, 0.75 and 0.50. From the numbers given there one deduces that the relative 
excess of H ± events is always quite large, especially at tan (3 = 30, and that the absolute 
statistics is huge, between O(10 4 — 10 5 ) signal events per year. In our opinion then, 
given the expected performances of the LHC detectors [p], 0, such signals could well be 
detectable soon after turning on the machine. Moreover, we stress that, as the total 
luminosity gets larger, the significance of the signal with respect to the total background 
will increase further. 

However, before drawing optimistic conclusions, one has to carefully consider first 
the kinematic properties of processes (1)— (4), as the LHC detectors will have a finite 
coverage (for example, in pseudorapidity r] and transverse momentum of the visible 
particles). Hence, we have plotted in Figs. 3 and 4, the differential spectra in the above 
variables for both signal and background processes. If one assumes that the phase space 
region that can be covered experimentally is approximately the one delimited by the 
requirements [77(6, 7", j) | < 3 and Pr(fe,r,j) > 20 GeV Jj], Q], then one can easily verify 
that most signal events are contained in the detectable region (we have checked that 
also in the case M4 = 60 GeV one gets distributions similar to those of Fig. 3 and 4) . 
This is true for background events as well, process (4) being possibly the only exception 
(in the pseudorapidity spectra, upper left corner in Fig. 4a and b). In the end then, 
one should expect that usual selection criteria will not alter the conclusions that were 

5 We believe that hadronic r decays can be easily distinguished from quark and gluon jets. 
6 In first approximation, we assume that the rejection factor for misidentification of light quarks 
and gluons as 6-jets is large enough, that we can neglect the reducible QCD background here. 



previously extracted from the total rates of Tab. I. 

If the charged Higgs mass is above the tb threshold, if ± scalars could reveal them- 
selves via the production and decay mechanism (5). For a heavy H boson, we consider 
then the signature bbb]]]X, where h represents either a quark or the corresponding an- 
tiquark, and j a jet that does not show a displaced vertex. In this case, the final state 
that should be detected is much more complicated than that of an intermediate mass 
boson, as it is made up of six jets. Nevertheless, the complex resonant structure 
of process (5) brings some advantages to the tagging procedure. In fact, on the one 
hand, three vertex tags are now required (which introduce the suppression factor ef), 
and in a high hadronic multiplicity environment; in addition, complications arise from 
the combinatorics of the the jets. On the other hand, the kinematics of the jets in the 
final state is highly constrained, since: (i) one of three possible jj combinations must 
reproduce the W ± mass (i.e., My ps M w ±); (ii) one of the nine possible b]] combina- 
tions must reproduce at the same time the t and W masses (i.e., M^j ps m t , with 
Mjj ps M w ±). In this respect, note that in the decay chain t — * bW ± — * 6jj we do 
not need to consider intermediate H^ 1 contributions, as in the heavy M#± range the 
t — > bH ± is forbidden. 

The total cross section for process (5) is displayed in Fig. 5, as a function of the 
A mass in the range 140 GeV ~ Ma ~ 480 GeV (for both tan/3 = 1.5 and 30). Both 
the t — > bW^ and — > jj branching ratios are included. For comparison, we also 
reproduce from Fig. 1 the rates for process (1). A first feature that is worth noticing 
in Fig. 5 is that again the rates of process (5) for tan (3 = 30 are larger than those for 
tan/3 = 1.5. In practice, there are two opposite effects which take place: on the one 
hand, at small tan /3, the 5.R(if ± — > tb) is larger whereas, on the other hand, at large 
tan /3, the contribution from graph 2 in Fig. la is enhanced by couplings proportional 
to tan/3 itselfQ. Of the two, it is the second that dominates over most of the M H ± 
range. Furthermore, at large tan/3's, processes (1) and (5) yield rates of the same order 
for Ma ~ 200 GeV, such that in this case they can be contemporaneously exploited 
in searching for H ± signals. This is not true at small tan (3. In general, it should 
be noted that in Fig. 5 one is dealing with total rates that are more than two orders 
of magnitude smaller that for the case of an intermediate mass if ± boson in the tv t 
channel. Nevertheless, as a starting point one can rely on O(10 — 100) signal events 
produced via — > tb per year (for the same luminosity as above). Clearly, in this 
case the need for a high value of is crucial. For example, a vertex tagging efficiency 
of 50% reduces the production rates by a factor of 8. 

In Fig. 6 we plot the spectrum in the sum of the invariant mass of all possible 
bb]] combinations entering in process (5)0. The values of Ma considered here as a 
reference are 200, 300 400 and 500 GeV. Since at least one of the bb]] systems is made 
up by the decay products of the charged Higgs boson, a peak should possibly appear 
in the distributions (at M H ± = 215,311,408 and 506 GeV, respectively), on top of the 
combinatorial background. Indeed, the resonant peaks are quite sharp (the H ± widths 

7 Also note that at the same time, since M H ± > m t + mj,, the first graph in Fig. la is no longer 
resonant. 

8 Note that in our NWA approach we automatically obtain M^j = m t for the right three jet 
combination. In the case of the W ± decay we have adopted the 'conservative' requirement — 
M^i^lSGeV. 



are approximately, for the above values of mass: 0.84, 3.03, 5.20 and 7.20 GeV at both 
tan /3 values), and clearly visible. However, the total number of events in the region, say, 
\M m -M H ±\ < 25 GeV is, for tan (3 = 30: 25(10)[3], 9(4)[1], 4(2) [0.4] and 2(1)[0.2] (for 
the four masses above), for 10 fb _1 per year of luminosity, assuming e b = 1(0.75)[0.5]. 
Therefore, at least for not too heavy H ±J s and high 6-tagging performances, one could 
possibly look for H ± signals in the tb channel at large tan /3's. Certainly, if the high 
luminosity option / Cdt = 100 fb _1 can be achieved at the LHC, things would be very 
optimistic, as in a few years of running, even the very heavy mass region could be 
scanned. At small values of tan/3 one has to consider rates that are typically smaller 
by one order of magnitude, rendering Higgs detection much more difficult. 

Finally, in Fig. 7, we show the differential spectrum in transverse momentum of 
the various bty] combinations that can be reconstructed from process (5), for the same 
values of Ma as above. We notice that from the figure it is clear that the pt spec- 
trum of the Higgs decay products is significantly hard (because of the large mass of 
the scalar), a feature that could well help in disentangling heavy H ± boson signals, 
especially considering that the ordinary QCD background in six-jet events has quite a 
soft transverse momentum distribution. 

4. Conclusions 

In this paper we have studied, within the AiSSAi, a new production mechanism of 
charged Higgs bosons H ± at LHC energies, via 6-quark initiated interactions. Two 
possible Higgs signatures have been considered, depending on whether the mass M H ± 
is below or above m t + m b ps 180 GeV. 

In the first case, by exploiting one vertex tag on one 6-jet in the final state, the 
Higgs decay channel H ± — > tv t should be identifiable as a clear excess in the number 
of r events with respect to the rates predicted by the non-SUSY backgrounds, provided 
that M A S 130 GeV (i.e., M H ± S 144 GeV), both at large and small values of tan/3. 
Also, the absolute number of signal events is statistically very large. In this mass range, 
a careful treatment of various background sources has been performed. 

In the second case, the channel H ± — > tb — > bb]] (via hadronic decays of the top) 
has been considered. The signature arising from heavy H ± decays is complicated (a 
six-jet final state involving three 6-jets) and only a small number of events is expected. 
Nonetheless, the resonant behaviour of the t, if ± and W ± decay products should allow 
one to eliminate the ordinary QCD background in light quark and gluon jets, although 
(in the heavy mass range) a signal-to-background analysis has not been performed. 
In general, if high 6-tagging performances can be achieved and/or the high luminosity 
option becomes available at the LHC, scalars with masses up to 500 GeV could 
well be searched for, as the combinatorial background does not spoil the form of the 
Higgs peaks. This is however true only for large values of tan/3, since for low tan/3's 
the event rates are smaller by one order of magnitude. 

The range between Ma ~ 130 — 140 GeV and up to the opening of the tb decay 
threshold is extremely difficult to cover, as rates in the rv T channel drastically decrease 
well below the background rates and at the same time the off-shell tb channel has a 
very small statistics. 



Finally, we stress that, before drawing any firm conclusion from our results, one 
should include a realistic simulation of the expected performances of the LHC detectors 
and that, in the heavy mass range Mh± > m t + nib, a detailed background study 
(including all the hadronisation effects in a six-jet final state, an analysis which was 
beyond our capabilities) should be performed. Nevertheless, we believe that the matter 
presented here would deserve experimental attention when proceeding to the various 
simulations of the M.SSM. Higgs phenomenology at the CERN hadron collider. 
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Appendix 

In this additional section we write down in analytic form the matrix element for the 
signal processes. As an example, we reproduce that of the reaction bU — > bDH + — ► 
bDr + v T © C.C.. However, by replacing r — > t and v T — > b one can easily obtain 
that for top-bottom (on-shell) production. In fact, we have used a Narrow Width 
Approximation for the top when this is produced from the H ± splitting, by computing 
the exact amplitude squared for bU — > bDH + — > bDtb © C . C . and by interfacing this 
with a SUBROUTINE implementing the top decay formulae into bW ± pairs as given in 



Ref. 0. 



The matrix element squared (summed/averaged over the final/initial spins and 
colours) for the process bU — > bDH + — > bDr + v T reads as: 

\M\ 2 = \M \ 2 \P H \ 2 \Pw\ 2 (\M t \ 2 + IM^l 2 + 2\M* t M^\), 

with 

\M \ 2 = (gV2M^) 2 [ PT - P , T (mlt a n 2 p + ml T Cot 2 p)-2m 2 T ml T } 
P H = (p 2 H -M 2 H + iM H T H )- 1 
Pw = (Pw ~M 2 V + iMwTw)- 1 
\M t \ 2 = 2\P t \ 2 p bM - Pu x 
x [2m 2 b (m 2 t + p btOUt ■ p t tan 2 (3) p t ■ p D + {m\ cot 2 (3 - m\p 2 t tan 2 (3) p bjOUt ■ p D ) 
\M^\ 2 = mlsec 2 (3(2p H ■ p D p H ■ Pu - p 2 H pu ■ Pd)* 
x[m 2 b (\P HoM \ 2 - \P Ao \ 2 ) + Pb,out ■ Pb,in (\Ph m? + \ p A \ 2 )] 
2l.Mf.M0l = m 2 sec/?x 
x [[ft[P t *(P HoM - P Ao )}m 2 cot (3 - ft[P;(P HoM + P Ao )]ml tan/5] x 
x (p b ,out -PdPh -Pu+ Pb.out -PuPh-Pd- Pb,out -PhPu ■ Pd) + 



+[U[P*(P Ho , ho - P Ao )\m 2 b tan (3 + U[P*(P HoM + P Ao )](m 2 t cot/3 + 2 Pbtin ■ p b)(mt tan/5)] x 

X (p 6 ,j n -PDPH-PU+ Pb,in ' PU PH ' PD ~ Pb,in ' PH PU ' Pd) + 

2$l[P t *(P HoM + P Ao )}t!mf3x 

X (Pb,in ■ PU Pb,out ■ PHPU ■ PD - Pb,in • PU Pb,out • PD PU " PH + 
+Pb,in -PdPV Pb,out PVPH- Pb,in ' PH PU " Pb,out PU ' Pd) ~ 2 ^^\aP bM Pb,outPuPD X 

x [S(P;P A) )m 6 2 tan /3 - Z(P t * P HoM )m 2 t cot/5 
-9f[P t *(-Pfl- ,ho + Pao)] tan/3(p 6>in • p h)OUt - p hM ■ pu + PbjOUt ■ pu)\. 
Here we have set m v = m D = 0, p t = p bjin + pw = Pb,out + Ph, and 

p t = (Pt ~ m t + im^t)' 1 
PhoM = cosasin(/3 - a)(pj - m 2 Ho + iM Ho r Ho )~ 1 + 
+ sinacos^ - a)(pj - m£ + i/l4 o r ho ) _1 
P A) = sin/5(4 - + iMa^aX 1 

PH = Pt+ Pu t P<f> = Pb,in - Pb,out Pw = Pu - Pd, 

where pb,in, Pu, Pb,out,PD,Pr and p Vr are the external momenta (incoming/outgoing in the 
initial/final state). Note that the symbols K and S refer to the real and imaginary part 
of a complex number, respectively, and that e^ u \ a is the Levi-Civita tensor (eoi23 = 1)- 
Finally, by exchanging p bjin <-> —p btOUt and e^ v \ a <-> — e^xa one can obtain the amplitude 
squared for bU fusion, and by relabelling U <-> P those for the 6P- and 6P-initiated 
reactions. 
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Table Captions 

[I] Number of events for signal (1) and total background (2)-(4), at the LHC with 
CM energy ^7s pp = 14 TeV, for M A = 60(80)[100]{120} GeV and tan/3 = 1.5 and 
30, assuming the yearly luminosity J Cdt = 10 fb _1 and the 6-tagging efficiencies 
Eb = 1.0,0.75 and 0.5 (first, second and third row, respectively), together with 
the relative excess of signal events respect to the total background. 



Figure Captions 



[1] Feynman diagrams contributing at lowest order to the subprocess bU — > bDr + v T ® 
C.C.: (a) via H ± resonant graphs; (b) via resonant graphs and through 
EW only; (c) via W ± resonant graphs and through QCD interactions; and to the 
subprocess UD — > bbT + u T @ C.C.: (d) via W ± resonant graphs and through QCD 
interactions. Here, U and D represent any of the light quarks u, d, s and c. The 
PostScript version of the Feynman graphs has been produced using MadGraph 

[2] Cross section of the four processes (l)-(4), at the LHC for \fs pp = 14 TeV, in 
the region 60 GeV ~ Ma ~ 480 GeV and for two different values of tan /3 (in the 
case of processes (1) and (2)). In the top left hand plot we enlarge the region 
60 GeV ~ Ma ~ 150 GeV. Solid line: process (1) with tan/3 = 1.5. Dashed line 
(large spacing): process (1) with tan/3 = 30. Dotted line (large spacing): process 
(2) with tan/3 = 1.5. Dot-dashed line: process (2) with tan/3 = 30. Dashed line 
(small spacing): process (3). Dotted line (small spacing): process (4). The cut 
pfinai > Q e y b een applied to all particles in the final states, except for 
neutrinos. The structure function set MRSA has been used. 

[3] Differential distribution in transverse momentum of the 6-quark (a) and of the 
r-lepton (b) for the four processes (l)-(4), at the LHC for \fs pp = 14 TeV, for 
the following selection of masses Ma = 80, 100, 120 GeV and for two different 
values of tan/3 (in the case of processes (1) and (2)). Solid line: process (1) with 
tan/3 = 1.5. Dashed line (large spacing): process (1) with tan/3 = 30. Dotted 
line (large spacing): process (2) with tan/3 = 1.5. Dot-dashed line: process (2) 
with tan/3 = 30. Dashed line (small spacing): process (3). Dotted line (small 
spacing): process (4). Please note that the two curves corresponding to process 
(2) are practically indistinguishable. The cut p^ nal > 10 GeV has been applied to 
all particles in the final states, except for neutrinos. The structure function set 
MRSA has been used. 

[4] Differential distribution in pseudorapidity of the 6-quark (a) and of the r-lepton 
(b) for the four processes (l)-(4), at the LHC for \fs pp = 14 TeV, for the following 
selection of masses Ma = 80, 100, 120 GeV and for two different values of tan (3 (in 
the case of processes (1) and (2)). Solid line: process (1) with tan /3 = 1.5. Dashed 
line (large spacing): process (1) with tan/3 = 30. Dotted line (large spacing): 
process (2) with tan/3 = 1.5. Dot-dashed line: process (2) with tan/3 = 30. 
Dashed line (small spacing): process (3). Dotted line (small spacing): process 
(4). Please note that the two curves corresponding to process (2) are practically 
indistinguishable. The cut p^ al > 10 GeV has been applied to all particles in 
the final states, except for neutrinos. The structure function set MRSA has been 
used. 

[5] Cross section of processes (1) and (5), at the LHC for \fs pp = 14 TeV, in the 
region 140 GeV ~ Ma ~ 480 GeV and for two different values of tan/3. Solid 
line: process (1) with tan/3 = 1.5. Dashed line (large spacing): process (1) with 
tan/3 = 30. Dashed line (small spacing): process (5) with tan/3 = 1.5. Dotted 



line: process (5) with tan/3 = 30. The cut p| nal > 10 GeV has been applied to 
all particles in the final states, including the top decay products. The structure 
function set MRSA has been used. 

[6] Differential distribution in invariant mass of the bb]] systems for process (5), 
at the LHC for \^s pp = 14 TeV, for the following selection of masses Ma = 
200,300,400,500 GeV, for tan/3 = 1.5 (solid line) and 30 (dashed line). The cut 
pfinai > GeV has been applied to all particles in the final states, including the 
top decay products. Bins are 5 GeV wide. The structure function set MRSA has 
been used. 

[7] Differential distribution in transverse momentum of the bbjj systems for process 
(5), at the LHC for \/s pp = 14 TeV, for the following selection of masses M A = 
200,300,400,500 GeV, for tan/3 = 1.5 (solid line) and 30 (dashed line). The cut 
pfinai > GeV has been applied to all particles in the final states, including the 
top decay products. The structure function set MRSA has been used. 
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